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Abstract 

In this note, we propose an original extended formulation of Euler-Korteweg systems based on a generalization of 
the quantum Bohm potential identity. This new formulation allows to propose a useful construction of a numerical 
scheme with entropy stability property under a hyperbolic CFL condition. We also comment the use of the identity 
for compressible Navier-Stokes equations with degenerate viscosities. 

Resume 

Generalisation de I’identite de Bohm quantique et nouvelles formulations augmentees pour equations 
d’Euler-Korteweg. Dans cette note, on propose une formulation augmentee originate des systemes d’Euler- 
Korteweg basee sur une generalisation de I’identite dite du potentiel de Bohm quantique . La motivation principale 
est la construction de schemas avec stabilite entropique sous condition CFL hyperbolique du systeme d’Euler- 
Korteweg. On presente egalement quelques commentaires concernant I’utilisation de cette identite dans le cadre 
des equations de Navier-Stokes avec viscosites degenerees. 


Version frangaise abregee 

Dans cette note, nous introduisons de nouvelles formulations augmentees (sous forme conservative) 
du systeme d’Euler-Korteweg (l)-(3) en plusieurs dimensions d’espace. Selon le choix des coefficients de 
capillarite, ce type de systeme intervient dans la modelisation des melanges de type liquide-vapeur, des 
super-fluides ou de I’hydrodynamique quantique par exemple. Nous nous interessons ici a de nouvelles 
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formulations permettant de construire un schema numerique a stabilite entropique. Nous etendons ainsi 
les travaux recents de deux des auteurs qui consideraient le cas uni-dimensionnel. Nous traitons les 
termes dispersifs implicitement et donnons un resultat de stabilite entropique des schemas d’ordre 1 sous 
condition CFL hyperbolique. A titre d’illustration, nous presentons des resultats numeriques pour des 
films minces avec tension de surface modelises par les equations de Saint-Venant. 


1. Introduction 

In this paper, we introduce (new) extended formulations of the so-called Euler Korteweg system which, 
in several space dimensions, reads 

9tp-|-div(pu) = 0, (1) 

dt{gu) + div(pu 0 u) -I- Vp{g) = div(K), (2) 

where g denotes the flnid density, u the fluid velocity, p{g) the fluid pressure and K the Korteweg stress 
tensor defined as 

K= (^gdiv{Kig)Wg) + ^{K{g)-gK\g))\Wg\^'^lM.-Kig)Wg(S)Wg. (3) 

with K(g) the capillary coefficient. These models comprise liquid-vapor mixture (for instance highly pres¬ 
surized and hot water in nuclear reactors cooling system) [7], superfluids (Helium near absolute zero) [6] 
or even regular fluids at sufhciently small scales (think of ripples on shallow waters) [8]. In quantum hy¬ 
drodynamic, the capillary coefficient is chosen so that gK{g) = constant: in this case; the Euler-Korteweg 
equations correspond to the nonlinear Schrodinger equation after Madelung transform. In classical fluid 
mechanics, the capillary coefficient K{g) is chosen constant. 

The system (l)-(2) admits two additional conservations laws. One conservation law is satisfied by the 
fluid velocity 

5tu-I-u • V u + V(5£) = 0, (4) 

with S the potential energy and S£ its variational gradient 

£{g.Vg)=F^{g) + ]^K{g)\\/g\\ 5£ = Fl,{g)-]^K'{g)\V g\^ - K{g)^g. (5) 

The local existence of strong solution to (1), (4) is proved in [1]. For that purpose, the authors introduced 
an extended formulation by considering an additional velocity w = V(/?(p) with y/gp'{g) = yjK{g): 

dtu -F u • V u + V ^Eg(£») - = V (a(p)div(w)), dtw + V (u^ w) = -V (a(£i)div(u)), (6) 

where a{g) = \/K{g)g. This formulation is particularly adapted to the derivation of a priori estimates, 
the very first one being a conservation law on the total (kinematic-|-potential) energy : 

(||up-h£(p, Vp)) -fdiv (^u(^||up+£:(p,V£<)-fp(p))) = div(^E(p)(Vw u - Vu w)^ 

—div^(F(£») — pE'(£»))(div(w)u — div(u)w)^. (7) 
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with F'{q) = q(p'{q). 


In this note, we introduce a new extended formulation of (l)-(2) by considering the conservative variables 
gw instead of u, w. The key point is the generalization of the quantum potential Bohm identity. It 
allows to transform the Euler-Korteweg system into a hyperbolic system perturbed by a second order skew 
symmetric term. The main motivation is the construction of a numerical scheme which is easily proved 
“entropy” stable. Two of the authors performed this approach in a one dimensional setting and proved 
entropy stability under “capillary” Courant-Friedrichs-Lewy (denoted CFL in the sequel) condition [9]. 
Here, we extend this approach to the multi-dimensional setting. Moreover, in order to avoid restrictive 
CFL condition, we treat dispersive terms implicitely and prove entropy stability of hrst order schemes 
under a hyperbolic CFL condition. We present preliminary numerical results for thin films with surface 
tension modeled by the shallow water equations. 


2. Generalization of the qnantum Bohm identity and extended formulation 

Let us first present an extension of the quantum potential Bohm identity 

2ediv(Ayg/v/g) = div(^£»VVlog£»^ 

strongly used in quantum fluid mechanics. More precisely, we can prove after some algebraic calculations 
the following relation 

e 

gV{./l^)AiJ =div(F(p)VV¥>(p))-v((F(^?)-F'(e)e)A(p(e)) ( 8 ) 

0 

with y/gip'ig) = yjK{g)^ F'{q) = \/K{g) g. This relation is a non trivial extension of the quantum Bohm 
identity which corresponds to the case K{g) = c/g. Remark that the left-hand side of ( 8 ) corresponds to 
the capillarity term K{g)\\/g\'^/2. It suffices to observe that 

Q 

K{g)\Vg\^ = \V{j 
0 

and thus as observed in [ 2 ] the variational gradient of the potential energy may be written 

6S = F^ig)-{y^A{J ^/Kis)ds)). 

0 

Following now the strategy of [1], we introduce a “good” additional unknown, homogeneous to a velocity. 
We denote this additional velocity w = ^ip{g) with y/g(p'{g) = \/K{g). In order to write a suitable 
extended formulation of the Euler Korteweg model, we also define F{g) so that F'{g) = \/K{g)g. The 
Euler Korteweg system admits the extended formulation 

/ 

dtg + div(£iu) = 0 , 

< dt{gn) + div(pu(g) u) -b Vp{g) = div(F(e)Vw^) - V {{F{g) - gF'{g))div{w)) , ( 9 ) 

9((pw) -b div(pw (g) u) = -div(F(p)Vu^) -b V {{F{g) - gF'{g))div{u )), 
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The total energy is then transformed into a classical entropy of the first order part of (9) 

||luf+f(p, Vp) = I (||uf+ ||wf) +Fo(p) :=£(p,u,w), 

whereas the second order part is skew symmetric. As a consequence, the energy conservation law (7) is 
obtained through a similar computation than in the first order case. 

Remark. Note that our formulation may be coupled to the result recently obtained in [5] to write an 
augmented formulation to the following compressible Navier-Stokes system with drag and capillary terms 

dtQ + div(pu) = 0, (10) 

dt{Qn) + div(pu 0 u) — 2div(^(p)I?(u) — V(A(p)divu) + Vp{q) + rip|u|“u = div(K), (11) 

if \{q) = 2{p'{q)q — p{g)) and K(g) = c{g.'{g))'^/g. Such compatible system could be used to prove 
global existence of weak solutions to the compressible Navier-Stokes equations with degenerate viscosities 
without capillary and drag terms. It has been recently performed in [10] in the case g,{g) = pg and 
X{g) = 0 introducing the quantum capillary term namely with K{g) = c/g. 


3. Application: stable schemes under hyperbolic CFL condition 


In this section, we introduce a numerical scheme for (9). The numerical domain is a rectangle defined 
by 0 < X < and 0 < y < Lj,, which is divided into fV = x riy rectangular cells. For the sake 
of simplicity, we consider uniform grid with constant spatial steps ix and by. We focus on the spatial 
discretization of the second order terms: they are written as da{f{g)dj 3 u) with {a, j3) G {x,y} and f{g) = 
F{g), gF'{g), F{g) — gF'{g). For that purpose, we introduce the following finite difference operators: 




/ , , ^i.j+1/2 ^i.j—1/2 / 7-1- \ /T \ 

{d2u),j = - {dJu),j+i/2= {d2u)ij = 


by 


by 


^iJ—1 
^ ' 


( 12 ) 


As a result the differential operator T{g)u = div(F(p)Vu^)-I-V {{gF'{g) — F(p))div(u)) is approximated 
by 7)i(p) defined as 


Fh(.g)^i.j — 


di {gF'{g)df ui)+ d2 {F{g)diU2)^ . + di {{gF'{g) - F{g))d2U2). . 
di {F{g)d2Ui)^^^ + 3,2 {{.gF'{g) - F{g))diUi)+ ^2 {gF'{g)d+U2)^ . 
We discretize (9) as follows 


/ ^n+1 n 


Q. 




el 


bt 

- (eu)r 


6t 


+ di{Fl,)i,,+d2iFl2kj = 0, 

+ di{Fl,k, + d2{Fl2kj = r,(p”+^)wK 


J ’ 


(pw)""!"^ — (pw)"- 

+di{K,^kj+d2{FZ,2k, = -r,(p”+i)uK+\ 


(13) 


where {k = 1,2) are classical Rusanov fluxes evaluated at p”, u”, w”. More precisely, the 

convection part is treated explicitly whereas the capillary terms are treated implicitly. Remark that there 
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are no capillary terms in the mass conservation law so that the implicit step amounts to solve a linear 
sparse system and is easily proved entropy stable. As a consequence, one can prove, by using discrete 
duality properties of the discrete second order operators, the following entropy stability result. 


Theorem 3.1 Suppose (13) is completed with periodic boundary conditions. Assume the hyperbolic scheme 
(system (13) with F = 0^ is entropy stable then the fully hyperbolic/capillary scheme (13) is entropy stable: 


2 = 0 ^=0 


,n+l 


,n+l __.n+l\ 


2 = 0 j—0 






This means that the numerical scheme (13) is entropy stable under a classical hyperbolic Courant- 
Friedrichs-Lewy condition. As an application, we carried out a numerical simulation of a thin film falling 
down an inclined plane. A consistent shallow water model [3] is given by 


dt{hu) + div (/lu (g) u) + V(p(/i)) + 


dfh + div {hu) = 0, 
dx Bi= S{h,\i) + yV(A/i). 


(14) 

(15) 


with p{h) = gcosl6)h^/2 and S')/;,, u) = ghsin(6*)ei — ?>v\i/h and ei the first vector of the canonical 
base directed downstream. Here g = 9.8 is the gravity constant, p, z/, a are respectively the fluid density, 
kinematic viscosity and surface tension whereas 6 is the inclination of the plane. We picked the values 
found in [8] for a solution with 31% glycerin by weight: p = 1.07 x 10^ kg.m~^, v = 2.3 x 10“® and 

cr = 67 X 10“^ kg.s~‘^. The source term is treated implicitly: since the source term is only in the equation 
for u and is linear with respect to u, the implicit step remains linear. We first carry out a numerical 
simulation of the original experience in [8] but imposed periodic boundary conditions in both directions. 



Figure 1. Numerical simulation of a roll-wave in presence of surface tension. On the left: one dimensional roll-wave without 
transverse perturbations. On the right: a two-dimensional roll-wave 


In order to test the robusness of the scheme, we also carried various numerical experiments of a drop 
falling down a plane in order to deal with wet/dry fronts. For that purpose, we introduced a precusor 
film with a thickness of 1.0 x 10“^ mm. 
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Figure 2. Drop falling down an incline plane (6 = 60°) at time £ = 0 and £ = Is. The fluid density, kinematic viscosity and 
surface tension are respectively p = 1.0 x 10^ kg.m~^, v = 1.0 x 10“® rn^.s”^ and cr = 67 x 10“^ kg.s~‘^. 

We will deal the problems of considering physical boundary conditions, deriving higher order schemes 
and considering wet/dry fronts in a forthcoming paper [4]. This will be useful to compute instabilities in 
moving contact lines. 
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